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Abstract 

Wootters [Phys. Rev. Lett. 80, 2245 (1998)] has given an explicit formula 
for the entanglement of formation of two qubits in terms of what he calls 
the concurrence of the joint density operator. Wootters's concurrence is de- 
fined with the help of the superoperator that flips the spin of a qubit. We 
generalize the spin-flip superoperator to a "universal inverter," which acts on 
quantum systems of arbitrary dimension, and we introduce the corresponding 
generalized concurrence for joint pure states of D\ x D2 bipartite quantum 
systems. We call this generalized concurrence the I-concurrence to empha- 
size its relation to the universal inverter. The universal inverter, which is a 
positive, but not completely positive superoperator, is closely related to the 
completely positive universal-NOT superoperator, the quantum analogue of 
a classical NOT gate. We present a physical realization of the universal-NOT 
superoperator. 



I. INTRODUCTION 

Entanglement plays a central role in quantum information theory |l]||. Perhaps the most 
important measure of entanglement for bipartite systems is the entanglement of formation 
Hi!]]. For a bipartite pure state \^ AB ), the entanglement of formation is given by the entropy 
of the marginal density operators, pa and pB, of systems A and B. For a bipartite mixed 
state pab, the entanglement of formation is given by the minimum average marginal entropy 
of ensemble decompositions of pab- 

Hill and Wootters introduced another measure of entanglement, called the concur- 
rence, for pairs of qubits. The concurrence is defined with the help of a superoperator S2, 
whose action on a qubit density operator p = ^(1 + P ■ a) is to flip the spin of the qubit: 
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S 2 (p) = a y p*a y = 1 -(I-P-a). (1.1) 

Here p* is the complex conjugate (or transpose) of p relative to the eigenbasis of o z . The 
concurrence of a pure state \^ AB ) of two qubits is defined to be || 



C 2 (^Ab) = ]J{^ABS 2 ®S 2 (\^ AB }(^ AB \) m AB ) = \(^ABWy®^ AB }\ . (1.2) 

Wootters |J showed that the entanglement of formation of an arbitrary two-qubit mixed 
state p AB can be written in terms of the minimum average pure-state concurrence, where 
the minimum is taken over all ensemble decompositions of p AB , and he derived an explicit 
expression for this minimum average pure-state concurrence. Wootters called this minimum 
the concurrence of the mixed state. Uhlmann |7j introduced a generalization of Wootters's 
concurrence to higher dimensions, which we discuss further in Sec. |11 B, 



In this paper we generalize the notion of concurrence to pairs of quantum systems of 
arbitrary dimension, in a way different from Uhlmann's. We show in Sec. |TJ that if the 
concurrence is to be generated by a product superoperator, as in the expression (|1.2|) , then 
the only suitable superoperator to go into the tensor product is what we call the "universal 
inverter." For a D-dimensional quantum system, which we call a "qudit," we denote the 
universal inverter by So- The action of the universal inverter on a qudit state p is given by 

S D (p) = v D (I - p) , (1.3) 

where is a positive constant. Acting on a pure qudit state the universal inverter 
maps \ip) to a multiple of the maximally mixed state in the subspace orthogonal to \ip). The 
universal inverter has been used previously in studies of the separability of mixed states by 
Horodecki and Horodecki [§]. 

The corresponding generalized concurrence for a joint pure state \^ AB ) of a D\ x D 2 
system, in analogy to Eq. (|1.2p for qubits, is given by 



C(m AB ) = yj{V A B S Di ®S D2 (\^ ab )(^ab\) ^ab) = 2u Dl u D . 2 [l - tr(p A )} . (1.4) 

Thus, for pure states, this generalized concurrence is simply related to the purity of the 
marginal density operators. A sensible choice for the constant Up, consistent with the 
concurrence for qubits, is v B = 1. We call the generalized concurrence (|1.4Q the I-concurrence 



to emphasize its relation to the universal inverter and also to distinguish it from a generalized 
concurrence introduced by Uhlmann 

The universal inverter is a natural generalization to higher dimensions of the qubit spin 
flip. Only for D = 2, the spin flip, does the universal inverter map pure states to pure 
states. The universal inverter cannot be realized as a quantum dynamics, because though 
it is a positive superoperator, it is not completely positive. In Sec. |I D| we explore a one- 
parameter family of trace-preserving superoperators that are closely related to the universal 
inverter, and we show that the completely positive member of this family that is closest 
to the universal inverter is the universal-NOT superoperator P,pH|]. The universal-NOT is 



thus a physically realizable quantum analogue of the classical NOT gate. The action of the 
universal-NOT, denoted {?not, on a qudit state is given by 
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In Sec. |TT| we give two physical realizations of the universal-NOT, one in terms of the 
quantum information distributor introduced by Braunstein, Buzek, and Hillery |TT| and the 
other in terms of a measurement of the isotropic POVM followed by state inversion. 

The paper concludes in Sec. ^V] with a brief discussion that includes the natural extension 
of I-concurrence to mixed states. 



II. UNIVERSAL INVERTER 



In this section we first review, in Sec. [11 A| , Wootters's spin-flip operation for a qubit and 
how it leads to an entanglement measure called the concurrence for an arbitrary pure state of 
two qubits || . The main result of this paper is to generalize the spin flip to a superoperator 
that we call the universal inverter. The universal inverter is defined in all Hilbert-space 
dimensions, and it leads to a generalized concurrence for joint pure states of two quantum 



systems of arbitrary dimension. In Sec. [H B| we formulate the requirements for the universal 
inverter and explore some of its properties, in Sec. [II Q we show that these requirements 
pick out a unique universal inverter up to a constant multiple, and in Sec. [11 D| we consider 
trace-preserving superoperators that are closely related to the universal inverter. 

The formalism we use for superoperators has been used extensively in open-systems 
theory ||12|| . The particular notation we use can be found in Ref. |13| and is summarized 
briefly in Appendix [A], along with a description of several superoperators that play key roles 
in our discussion. In contrast to Ref. fl3|| , we use 0, instead of <g>, to denote the slot into 
which one inserts the operator on which a superoperator acts, reserving ® to denote tensor 
products between quantum systems. This superoperator formalism has been used to analyze 
entanglement in Ref. [|14[ . 

We refer to the two subsystems of a bipartite system as systems A and B. Where 
necessary for clarity, we use subscripts A, B, and AB to distinguish quantities belonging 
to the subsystems and to the joint system. To reduce notational clutter, however, we omit 
these subscripts on pure states, denoting pure states of a single system by a lower-case Greek 
letter, e.g., \ifi), and joint pure states of a bipartite system by an upper-case Greek letter, 
e.g., I*). 



A. Spin flip and qubit concurrence 

A spin flip for a single qubit is effected by the anti-unitary operator a y C = —Co y , where 
C denotes complex conjugation in the eigenbasis of a z . Acting on a state vector or an 
operator A, the anti-unitary complex conjugation operator gives C\ip) = or CA = A*C, 
where and A* denote complex conjugation of the state or operator in the eigenbasis of 
o~ z . An anti-unitary operator that satisfies G 2 = ±1, i.e., Q> = ±B, is called a conjugation; 
conjugations are ordinarily introduced in quantum mechanics to represent time reversal. 
Complex conjugation in some orthornormal basis is a conjugation because C 2 = 1, and spin 
flip is a conjugation because (o~ y Cy = C^a y = Ca y = —a y C. For a description of other 



properties and uses of anti-linear operators, see Ref. [15 
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Promoted to an operator on operators, the spin flip becomes an anti-linear superopera- 
tor a y C C<Ty, which acts on operators according to a y CACa y = a y A*a y . Since we are only 
interested in the operation of the spin flip on Hermitian operators, where complex conju- 
gation is equivalent to transposition, we can replace this anti-linear superoperator with the 
corresponding linear superoperator 

S 2 = cr y o y o % , (2.1) 

where T 2 denotes transposition in the eigenbasis of a z (see Appendix [A]). The subscript 2 
distinguishes the spin flip and transposition in two dimensions from the similar quantities 
for arbitrary dimensions that we introduce later in this section. 

The action of the spin-flip superoperator on an arbitrary qubit density operator, p = 
I (I + P • a), is to invert the Bloch vector P through the origin, as in Eq. ( |1 . 1|) . Since 
inversion commutes with rotations, representing unitary operators, we have immediately 
that S 2 commutes with all unitary operators U, i.e., £2 o U ift = U G)WoS 2 . 

For a quantum state p of a two-qubit system, the spin-flipped density operator, distin- 
guished by a tilde, is 

P = S 2 S 2 (p) = Oy® CTyP*(Ty (Jy . (2.2) 

Hill and Wootters || defined the concurrence of a two-qubit pure state, p = \^/) to be 



C 2 (tt) = Vtr(pp) = yJ(V S 2 <g>Sj(|tf)<#|) *) = K*|(T W a y \**)\ . (2.3) 
The joint pure state can be written in terms of a Schmidt decomposition, 

|*> = ai\ei) \fi) + a 2 \e 2 ) |/ 2 > , (2.4) 

where \ej) and \fj) are the orthonormal eigenvectors of the marginal density operators for 
the two qubits and a± and a 2 are the (positive) square roots of the corresponding eigenvalues. 
Since S 2 commutes with all unitary operators, the concurrence C 2 (^) is unchanged by local 
unitary transformations. This means that C 2 {^) is a function only of a\ and a 2 \ it is 
easy to verify that C 2 (ty) = 2a\a 2 . As noted by Wootters, the concurrence can serve as a 
measurement of entanglement: it is invariant under local unitary transformations, as any 
good measure of entanglement should be, and it varies smoothly from for pure product 
states to 1 for maximally entangled pure states. 

Wootters [|] went on to show that the concurrence can also be used to measure the 
entanglement of mixed states of two qubits. He showed that the entanglement of formation 
of an arbitrary two-qubit mixed state p can be written in terms of the minimum average 
pure-state concurrence of ensemble decompositions of p, and he derived an explicit expression 
for this minimum in terms of the eigenvalues of pp. Wootters called the minimum the 
concurrence of the mixed state. 

Uhlmann 0] based his generalization of concurrence on the fact that the spin flip is 
a conjugation, defining concurrence in arbitrary dimensions in terms of a conjugation 0. 
Following Uhlmann, we call his concurrence the 0- concurrence, and we discuss it further in 
the next subsection. 
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B. Universal inverter and I-concurrence 



Our goal in this paper is to generalize the spin-flip superoperator S 2 for a qubit to a 
superoperator Sd that acts on qudit states and generates a generalized concurrence for 
Di x D 2 bipartite quantum systems. The spin-flip superoperator has several important 
properties that we might wish its generalization to retain: 

1. S 2 maps Hermitian operators to Hermitian operators. 

2. S 2 commutes with all unitary operators. 

3. (fy\S2 <S> i!>2(| x I / )( l I / |)| x I / ) is nonnegative for all joint pure states |^) and goes to zero if 
and only if \^f) is a product state. 

4. ^2 is a positive superoperator; i.e., it maps positive operators to positive operators. 

5. S 2 is trace preserving. 

6. S 2 maps any pure state to the orthogonal pure state \ip ± ){?p~ L \- 

7. S 2 is derived from a conjugation 0, i.e., an antiunitary operator satisfying 2 = ±1. 

Property 1 guarantees that S 2 (g) S 2 maps Hermitian operators to Hermitian operators (see 
Appendix |B|) and thus that the quantity (^\S 2 ® S 2 (\^)(^\)\^) of property 3 is real. Prop- 
erty 2 ensures that C 2 (^) is unchanged by local unitary transformations, as an entanglement 
measure should be. Property 3 makes C 2 (^>) well defined, by ensuring that the quantity 
inside the square root is nonnegative, and it sets the zero so that pure product states, but 
no other pure states, have vanishing concurrence. 

In generalizing the spin flip to higher dimensions, we want the generalized concurrence 
of a pure state p = \^/)(^\ of a D 1 x D 2 bipartite system to be defined as for qubits, i.e., 



It is clear that the analogues of properties 1-3 are desirable properties of Sd, for the same 
reasons as for qubits, and it turns out that they are sufficient to pick out a unique superop- 
erator Sd up to a constant multiple. 

The upshot of this discussion is that we require Sd to have the following properties: 

1'. Sd maps Hermitian operators to Hermitian operators. 

2'. Sd commutes with all unitary operators. 

3'. ('$\Sd 1 ® ^AjG^X^DI^) is nonnegative for all joint pure states \^) and goes to zero 
if and only if is a product state. 

The only superoperator that has these three properties is 




(2.5) 




(2.6) 
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where I is the unit superoperator relative to the left-right action, X is the unit superoperator 
relative to the ordinary action, and vd is an arbitrary real constant. For the considerations 
in Sec. [11 PL we allow vd to have a dependence on D. For purposes of defining a generalized 



concurrence, however, vd should be independent of D; otherwise the generalized concurrence 
of joint pure state could be changed simply by adding extra, unused dimensions to one or 
both systems. 

We show that Sd is the only superoperator allowed by properties l'-3' in Sec. [11 Q For 



the remainder of this subsection, we show that Sd does satisfy properties l'-3', and we spell 
out some of its other properties and properties of the corresponding concurrence. Notice 
first that takes an operator A to 

S D (A) = u D [I{A)-l{A)} = u D [tr{A)I - A) , (2.7) 

from which it is clear that Sd satisfies properties 1' and 2'. If A is a density operator p, we 
get 

S D (P) = M I -P) • (2- 8 ) 

Since I — p is a positive operator for any p, we have immediately that Sd is a positive 
superoperator provided that vd is positive. The generalized concurrence is indifferent to a 
change in the sign of vd, so we are free to choose vd to be positive, which we do henceforth, 
thus making Sd positive. If v D = 1/(D — 1), Sd is trace preserving; this trace-preserving 
normalization is useful for the considerations of Sec. Ill D|, but we see below that Vd = 1 is 



a more reasonable normalization to use for the generalized concurrence C(\l/). Finally, Sd 
maps a pure state p = I^K^I to a positive multiple of the projector orthogonal to p: 

s D mm = v D (i-\ip)m . (2.9) 

It is this property that prompts us to call Sd the universal inverter. We call the corre- 
sponding generalized concurrence Q2.5| ) the I-concurrence to emphasize its connection with 
inversion. Other properties of Sd, which follow directly from the corresponding properties 
of I and X (see Appendix [A]), are that Sd is Hermitian relative to the ordinary action, 
i.e., Sd = Sd, and that it changes sign under sharping, i.e., Sf) = —Sd- The universal 
inverter has been used previously by Horodecki and Horodecki || to provide a criterion for 
the separability of mixed states. 

We now see that properties 4-6 of the qubit spin flip survive, in amended form, in its 
generalization: 

4'. Sd is a positive superoperator. 

5'. Sd is a positive multiple of a trace-preserving superoperator, i.e., Sd(I) = ud(D — 1)1. 

6'. Sd maps any pure state \ip)(ip\ to a positive multiple of the projector onto the subspace 
orthogonal to 

It is worth pointing out that if we added to properties l'-3' the additional requirement that 
Sd map each pure state to a multiple of some orthogonal state, then the superoperator of 
Eq. (|2.6|) would trivially be the only possibility for the universal inverter. 
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We still have to deal with property 3'. For that purpose we need the tensor-product 
superoperator 



S Dl ®Sd 2 = v Di vd 2 {1®I- J®I-I®X + . (2.10) 

Applied to an arbitrary joint density operator pab, this tensor-product superoperator gives 

S Dl ® Sd 2 (pab) = v Dl v D2 {I ® I - p A ® I - I ® p B + p AB ) ■ (2.11) 

Projecting back onto pab gives 

h{pabS Di ® S d . 2 ( P ab)) = u Dl u D2 [l - tr(p A ) - tr(p|) + tr(p^)] > . (2.12) 

The inequality here, which shows that the quantity in property 3' is nonnegative, is proved 
in Appendix [C], where it is also shown that the inequality is saturated if and only if pab = 
Pa ® Pb is a product state, with pa or ps a pure state. For a joint pure state pab, this 
establishes property 3'. 

It is useful to specialize Eq. ( |2.12| ) to a joint pure state in which case it becomes the 



square of the pure-state I-concurrence: 

c 2 m = (* s Dl ® s D2 mm *) = 2^ 1 ^ 2 [i - tr^)] . (2.13) 



Thus the I-concurrence measures the entanglement of a pure state in terms of the purity, 
t r (pi) = tr(p^), of the marginal density operators. A joint pure state has a Schmidt 
decomposition, 



|*> = X)a i |e i )®|/ i ), aj >0, (2.14) 
j 

in terms of which the squared I-concurrence becomes 

C 2 ^) = 2v Dl v D2 (l-Es 4 ) =4z/ Dl ^ 2 Es 2a fc- ( 2 - 15 ) 
V j ) j<k 

For defining a concurrence, one should choose the scaling factor vd to be independent of D— 
otherwise, as noted above, the pure state I-concurrence could be changed simply by adding 
extra, unused dimensions to one of the subsystems — and to be consistent with the qubit 
concurrence, one should choose ujj = 1. With this choice the pure-state I-concurrence runs 



from zero for product states to y2(M — 1)/M, where M = min(Di, D 2 ), for a maximally 
entangled state. 

Of the seven properties of the spin flip listed above, the first six survive, some in amended 
form, in the universal inverter. The seventh, that S2 is derived from a conjugation, is not 
a property of So, because except in two dimensions, a conjugation cannot commute with 
all unitaries and thus cannot serve as the basis for a measure of entanglement. Uhlmann's 
work on conjugations is valuable in that it generalizes to all conjugations the expression 
that Wootters gives for the minimum average pure-state concurrence of a bipartite density 
operator p in terms of the eigenvalues of pp. Our results show, however, that Uhlmann's 
O-concurrence 0, founded as it is on the use of conjugations, cannot serve as the basis for 
a general measure of entanglement. 
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There is another interesting form of the universal inverter, which makes a direct con- 
nection to the form (|2.2| ) of the spin flip. Choosing an orthonormal basis \ej), let T be 
the superoperator that transposes matrix representations in this basis, and let Va be the 
superoperator projector, relative to the left-right action, which projects onto the subspace 
of operators that are antisymmetric in this basis. We show in Appendix [A] that 

S D /v D = 2V A oT . (2.16) 

This form of the universal inverter has been given previously by Horodecki and Horodecki 
||. For qubits, if we use the eigenstates of a z as the chosen basis, then the antisymmetric 
operator subspace is spanned by the normalized operator a y j a/2, so the projector onto this 
subspace is Va = l cr y)( cr ?/|/2 = cr y Qa y /2. Thus in the two dimensions the universal inverter 
becomes S2 = V2&y &y which agrees with the spin flip if v-i = 1. 



C. Derivation of universal inverter 

We now show that the only superoperator that satisfies properties l'-3' of the preceding 
subsection is the universal inverter ( |2.6|) . As we proceed through the proof, we use Qd to 
denote the operator under consideration. 

As we show in Appendix ||, property 1' implies that Qd is left-right Hermitian, i.e., 
Qd = Qdi an d thus has an eigendecomposition 

Qd = J2^ a \r a )(T a \ = ^2fi a r a d)Tl , (2.17) 

a a 

where the \x a are real (left-right) eigenvalues and the operators r Q are the corresponding 
orthonormal eigenoperators. 
Property 2' implies that 

Qd = C/ f U o Q D o U © t/t = J2t*aU j T a U iPr^U , (2.18) 

a 

which means that Wr a U is an eigenoperator of Qd, with eigenvalue fi a , for any unitary 
operator U . This result can be restated as saying that the degenerate eigensubspaces of 
Qd are invariant under all unitary transformations. We show in Appendix |B| that the 
only operator subspaces that are invariant under all unitary transformations are the one- 
dimensional subspace spanned by the unit operator and the (D 2 — l)-dimensional subspace 
of tracefree operators. As a consequence, Qd must have the form 

Q D = I2 D 1/D + v d F . (2.19) 

Here X — I I is the unit superoperator relative to the ordinary action, T is the superop- 
erator that projects onto the subspace of tracefree operators when acting to the right (see 
Appendix 0), \id is the eigenvalue of Qd corresponding to the normalized eigenoperator 
I/yD, and Vd is the eigenvalue corresponding to all of the tracefree operators. Notice that 
Qd is Hermitian relative to the ordinary action, i.e., Qd = Qd- 

If we add I/yD to a complete, orthonormal set of tracefree operators, we obtain a 
complete, orthonormal set of operators, so the unit superoperator in the left-right sense is 
given by 
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I = I/D + T , (2.20) 

from which we get 

Q D = VD l + v D \ , (2.21) 

where 

Vd = {»d ~ v D )/D . (2.22) 

Now we impose property 3'. In doing so, it is sufficient to consider the requirements of 
property 3' in the case where the two subsystems have the the same dimension D. In this 
case the tensor-product superoperator takes the form 

Gd®Gd = Vd' i ®' i + VdVd(Z ® I + I ® X) + vjjl ® I • (2.23) 
Applying this superoperator to a joint density operator pab gives 

Gd ® Gd(pab) = V 2 dPab + Vdv D (pa ® I + I ® Pb) + v\>l ® I , (2.24) 
and projecting this back onto pab yields 

^[pabGd ® Gd(pab)) = Vd^ t (Pab) + VD^D[tr(p 2 A ) + tr(p|)] + v 2 D . (2.25) 
Specializing to a joint pure state we get 

^\Gd ® |)|*) = + + 2^z/ D tr(p^) = (r/ D ^ z/ D ) 2 ± 2^ D z/ D [l ± tr(p^)] . 

(2.26) 

If tjd^d > 0, the top sign in Eq. ( |2.26|) shows that the quantity in property 3' is strictly 
positive, unless r] D = u D = 0, a case of no interest. If rj D u D < 0, the bottom sign in Eq. (]2.26|) 
shows that the quantity is nonnegative and goes to zero if and only if rjjj = —vd and pa 
is pure, i.e., the joint pure state is a product state. Thus it turns out that the quantity in 
property 3' is nonnegative for all superoperators of the form ( |2.21| ), but the only way to set 
the zero properly is to choose r)r> = —vd, thus giving the universal inverter of Eq. (|2.6j ). The 
left-right eigenvalues of the universal inverter are u D and po = Drj D + v D = —{D — l)v D . 

D. Trace-preserving superoperators 



All superoperators of the form (|2.21|) are proportional to a trace-preserving superopera- 
tor, since 

GW) = Gd{I) = (Vd + Du D )I . (2.27) 
Requiring Gd to be trace preserving gives the condition 

VD = 1- Dv d (2.28) 
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[po = D — vd{D 2 — 1)], which allows us to eliminate one parameter and to write the 
trace-preserving version of Qo as 



g DT = {\- Du D )X+u D l. (2.29) 

Acting on an arbitrary input state p, this superoperator gives 

G DT (p) = (1 - Dv D )p + v D I . (2.30) 

It is instructive to investigate this one-parameter family of trace-preserving operators. 

We first ask which of the trace-preserving operators ( |2.29|) are completely positive. The 
condition that a superoperator be completely positive is that its left-right eigenvalues be 
nonnegative (see Appendix |A]). Thus the condition for the complete positivity of Qdt is 
that pd > and vd > 0, which is equivalent to 

< u D < . (2.31) 

When u D = 0, Q DT = X is the unit superoperator, and when v D = D/(D 2 — 1), 

° J r =^~-(DI-l) = g NO T (2.32) 



WT 



D 2 -l D 2 -l 



is the universal-NOT superoperator |p|,|T0|| . Notice that the universal-NOT is a multiple of J 7 , 
the superoperator whose right action projects onto the subspace of tracefree operators. Since 
the dynamics of a quantum system must be completely positive, the universal-NOT is the 
closest physical approximation to the universal inverter in the one-parameter family ( |2.30P ; 



it is the quantum analogue of the classical NOT gate. We present a realization of the 
universal-NOT in Sec. pTT[ 

Another interesting completely positive superoperator occurs for ujj — + 1): 

Qdt = — ^(1 + X) = i J + — = G AV . (2.33) 

This superoperator was used to generate operator expansions in Ref. , where it was shown 
that it is the unique trace-preserving superoperator that satisfies Q = = Q x = and 
commutes with all unitaries. In contrast, the universal inverter is the unique superoperator 
that satisfies Q = = Q x = —Q* and commutes with all unitaries. 

As shown in Ref. fl4] , the superoperator Qay is the trace-preserving version of the su- 



peroperator that describes projection onto a random pure state, 

Gav = dJ^-\^)^\Q^)^\, (2.34) 

where dV is the unitarily invariant integration measure on projective Hilbert space and V is 
the corresponding total volume. Projection onto a random pure state is the measurement 



that results in the optimal estimation of the state of the qudit ||16|| . This estimated state is 
given by the density operator 

£av(p) = + P) ■ (2.35) 
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The superoperator Qav returns in Sec. [Hj as an ingredient in one of the physical realizations 
of the universal-NOT. 

We now consider which of the trace-preserving operators ( |2.29| ) are positive. Letting pj 
be the eigenvalues of the input density operator p, one sees that the eigenvalues of Qdt(p) 
[Eq. [2.30D 1 are (1 — Dv^Vj + VD- The condition that Qdt be positive is that these eigenvalues 
be nonnegative for all input eigenvalues pj, which is equivalent to 

< u D < ^—j . (2.36) 

When vo = 1/(D — 1), Qdt becomes the trace-preserving version of the universal inverter, 

The positive superoperators are convex combinations of X and Sdt'- 

Qdt = [1 - M D ~ + MD - l)S DT • (2.38) 
Notice that the universal-NOT can be written as 

£not = -^{Sdt + Qav) ■ (2.39) 



III. PHYSICAL REALIZATIONS OF THE UNIVERSAL-NOT 

In this section we give two physical realization of the universal-NOT superoperator (?not 
of Eq. Q2.32j ), the first in terms of the quantum information distributor introduced by Braun- 



stein, Buzek, and Hillery [TT[ and the second in terms of a measurement of the isotropic 



POVM followed by state inversion. 

For the first, consider a qudit in a pure state p = 1-0} (^l- As shown in Sec. [D], the ideal 
inversion of this state is given by 

Sdt(p) = -^y(/ - p) = / , (3.1) 

where Sdt is the trace preserving version of the universal inverter [see Eq. ( |2.37| )]. The 
inverted state p L is the maximally mixed state in the (D — l)-dimensional subspace orthog- 
onal to the input state p = \ip)(ip\. Notice that by construction, tr(pp ± ) = for pure input 
states. 

As shown in Sec. [II D| , the trace-preserving universal inverter 5dt is a positive, but not 



completely positive superoperator and as such cannot be realized physically. In the one- 
parameter family of trace-preserving inverters considered in Sec. [II D , the universal-NOT 



superoperator £not of Eq. ( |2.32| ) is the closest completely positive superoperator to the 
universal inverter. We denote the physically possible inversion of the state p obtained using 
the universal-NOT as 

p not = Wp) = d^ {di ~ p) ■ (3 ' 2) 



n 



In order to realize the universal-NOT, we couple the qudit to be inverted, denoted by A, 
to the quantum information distributor (QID) introduced in Ref. The QID is composed 
of two ancilla qudits, B and C, each of which has the same dimension D as qudit A. To 
describe the universal inverter, we introduce several operators and states for qudits. 

First we need the conjugate "position" and "momentum" operators, x and p. The 
eigenvectors of x are denoted by \xk), 

x\x k ) = x k \x k ) , (3.3) 



with the eigenvalues given by x k = k^lixjD; analogously, the eigenstates of p are denoted 
by \Pk), 

p\Pk) = Pk\Pk) , (3.4) 



with the eigenvalues given by p k = k^2ir/D. We use units such that the two operators 
are dimensionless. The two sets of eigenvectors, {!#*;)} and {\pk)}, form bases in the qudit 
Hilbert space and are related by a discrete Fourier transform, 



\xk) = ^=T,e-^ kl / D \ Pl ) , (3.5) 

v v 1=0 

W = 4E^ /D W- (3-6) 



D k=o 

The translation (shift) operators, defined by 

R x (n) = e - lx "P , R p {m) = e iPmX , (3.7) 

cyclically permute the basis vectors according to 

Rx(n)\x k ) = |x( fc+n ) mod£ )) , (3.8) 
R p (m)\pi) = \p(i +m )modD) , (3.9) 

where the sums of indices are taken modulo D. 

An orthonormal basis of D 2 two-qudit maximally entangled states |5 mn ) is given by 

\E mn ) = ^= Dy £e 2mmk l D \x k ) ® |x (fc+n)modD ) , (3.10) 

V L) k=0 

where m, n = 0, . . . , D — 1. Using Eq. ( |3.5| ), we can rewrite the states |H mn ) in the joint 
momentum basis: 

i D-l 

\E mn ) = -= e- 27Tml/D \p (m -i)modD) ® \pi) ■ (3.11) 

V L) t=0 

The state |S 00 ) can be written as 

1 , 

|2 00 ) = -7= \ X k) ® \ X k) = -7^ \P-lmodD) <S> \pi) ■ (3.12) 

V±> k=n vD z=n 



2 D-i y D ~ l 
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It is interesting to note that the whole set of D 2 maximally entangled states |S mn ) can be 
generated from |S o) by the action of local unitary operations (shifts): 

|S m „) = R p (m) <g> R x (n)\E 00 ) • (3.13) 

Now we are ready to describe the QID. The ancilla qudits, B and C, are initially 
prepared in the state 

|$) BC = 6|5 o)bc + &\x ) B ® |p )c • (3.14) 

The phase freedom in |$) BC can be used to make £1 real and nonnegative, but then £ 2 is in 
general complex. We do not use the freedom to make £ x nonnegative, thereby retaining for 
use below the ability to multiply both £1 and £ 2 by —1. 
Normalization of \&bc) imposes the constraint 

1 = £ + I6I 2 + ^ = £ 2 + ^ + b 2 + ^ , (3.15) 

where £2 = + Solving for £1, we get 



a / D 2 - 1 

a—Tj + f-^-a 2 -^- ( 3 - 16 ) 

We discard the other solution of the quadratic equation, because it can be converted to this 
solution by multiplying both £1 and £2 by —1. Since £1 is real, we must have 

^^a 2 + 6 2 <l, (3.17) 

which means that £2 lies on or within an ellipse that has principal radius D j \/D 2 — 1 > 1 
along the real axis and principal radius 1 along the imaginary axis. Therefore, we conclude 
that 

0<|6| 2 <^r. (3.18) 

It is easy to see that the minimum value of £1 occurs when £2 = D / \JD 2 — 1, this minimum 
value being £ x = —\j\jD 2 — 1. It is also easy to see that the maximum value of £1 occurs 
when £2 is real; the maximum occurs at £ 2 = — 1/V-D 2 — 1 and is given by £1 = D/y/D 2 — 1. 
The upshot is that £1 is bounded by 

1 D 

<&<-__. (3.19) 



The negative values of £ x are unimportant, because they can be converted to positive values 
by multiplying both £1 and £2 by —1. What is important is that |£i| 2 has the same range of 
possible values as |£ 2 | 2 - 

We now allow qudit A to interact with the two ancilla qudits, the resulting dynamics 
described by the unitary operator 
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Uabc = exp[-i(x c - x b )pa] exp[-ix A (pB + Pc)] (3.20) 



(for more details, see Ref. [ill). For an initial pure state of qudit A, the joint state after 
the interaction is 

UabMa ® \&)bc = ZMa ® \^oo)bc + 6 Mb ® |2 00 )ac • (3.21) 
The output states of the individual qudits after tracing out the other two qudits are 

pr ) =(e i 2 + ^ Z ^)p + % ! /, (3.22) 

p r ) =(N 2 + ^^)p + |/, (3.23) 



(out) _ 6(6 + 6>) r i £i + I6I 2 T 



where p is an arbitrary initial state of qudit A and p T is its transpose. Taking into account 
the constraint ( ft.15 ), we can rewrite the output states of qudits A and B as 

P 5 ut) = (l-|6| 2 )p+l6| 2 /AD, (3.25) 

p£ ut) = (i-6 2 )p + 6 2 '/£>- (3-26) 

As far as qudit A is concerned, the QID acts like the superoperator Q DT of Eqs. Q2.29Q and 
( p.30| ) with Dv D = |£ 2 | 2 - As far as qudit B is concerned, the QID first swaps the states of 
A and B and then acts like Qdt with Dv^ = £f. 

Rewriting the output state of qudit A in terms of the ideal inverted state p x = (I — 
p)/(D - 1), we get 

p%*> = (|6| 2 - 1)(D - 1)/ + [D — |6| 2 P - 1)]I/D . (3.27) 
To make p^^ as close as possible to /r 1 , we need to maximize |^2| 2 ; i.e., we need to choose 

Dv D = I6I 2 = j^r~[ > ( 3 - 28 ) 

thus making the action of the QID on qudit A the same as the action of the universal-NOT 
given in Eq. (|3.2|) . Notice that the QID gives the superoperator Qav of Eq. (|2.33|) when 
Dv D = \Z 2 \ 2 = D/{D + 1). 

When |£ 2 | 2 has its maximum value, £f — 1/ (D 2 — 1), so the output state (|3.26|) of qudit 



B becomes 

Notice that in the limit of large D, we have |£ 2 | ~~ *■ 1 an d 6 ~~ > 0- The output state of qudit 
B reduces to the input state of qudit A, and the output states of A and C reduce to the 
maximally mixed state I/D. All this is a consequence of the fact that the initial state of 
qudits B and C limits to \Q)bc - > \ x o) b ® \Po)c, an d the QID swaps the states of A and B: 
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Uabc\^)a ® |S o)sc = W)b ® |5 00 )ac • (3.30) 

Our second realization of the universal-NOT starts with a measurement of the isotropic 
POVM 

dE(m = D^-\i;)^\, (3.31) 

where 

r r dV 

J dE{m = Dj = (3.32) 

We assume that the measurement projects the system onto the measured state, so the 
operation that describes a measurement whose result is the state \ip) is 

dA(m = D^-\^\Q\^)(rlf\. (3.33) 

Knowing that the system is in the state we can invert the state. The operation that 
describes the measurent followed by inversion is Sdt ° dA(\ip)), where S^t is the trace- 
preserving version of the universal inverter. If we now throw away the result of the mea- 
surement of the isotropic POVM, the resulting trace-preserving operation is 

J S DT o dAQi/))) = S DT o £?av , (3.34) 

where Qay is the superoperator that describes projection onto a random pure state [see 
Eq. ( gap ]. 



Using the forms (|2.37|) and (|2.33|) , we can write the overall operation as 



S DT o Gay = jJ—ril - 2) o (I + 1) = -J—r(DI - 1) = £ NOT , (3.35) 
D A — 1 D A — 1 

where we use the fact that I o I = DI. This demonstrates that the universal-NOT results 
from a measurement of the isotropic POVM followed by state inversion. 



IV. CONCLUSION 

The concurrence introduced by Hill and Wootters and by Wootters || provides a 
good measure of the entanglement of any state of two qubits, pure or mixed. The Hill- 
Wootters concurrence is generated with the help of the superoperator that flips the spin of a 
qubit. In this paper we have identified the crucial properties of the spin-flip superoperator, 
which allow it to generate a good entanglement measure for pure states of two qubits. 
By generalizing these properties to systems of arbitrary dimension, we have singled out a 
unique superoperator, which we call the universal inverter. In the same way that the spin flip 
generates a concurrence for pairs of qubits, the universal inverter generates a concurrence, 
which we call the I-concurrence, for joint pure states of pairs of quantum systems of arbitrary 
dimension. This pure-state I-concurrence measures entanglement in terms of the purity of 
the marginal density operators of the joint pure state. 
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It is natural to define the I-concurrence of mixed states of D\ x D2 quantum systems as 
the minimum average I-concurrence of ensemble decompositions of the joint density operator. 
Property 3' of the I-concurrence — that the I-concurrence of a pure state is zero if and only 
if \ty) is a product state — implies immediately that the mixed-stated concurrence just defined 
is zero if and only if the mixed state is separable. We are investigating further properties 
of this mixed-state I-concurrence and how it is related to other measures of mixed-state 
entanglement. 

The universal inverter turns out to be the ideal inverter of pure states, since it takes 
a pure state to the maximally mixed state in the subspace orthogonal to the pure state. 
Because the universal inverter is a positive, but not completely superoperator, it cannot 
be realized as the dynamics of a quantum system coupled to an ancilla. We have shown 
that among a one-parameter family of inverting superoperators, the completely positive 
superoperator that comes closest to achieving an ideal state inversion is a superoperator 
called the universal-NOT, and we have presented a physical realization of the universal- 
NOT. 
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APPENDIX A: SUPEROPERATOR FORMALISM AND SPECIAL 

SUPEROPERATORS 

The formalism we use for superoperators has been used extensively in open-systems 
theory ||I2 |. In this Appendix, we summarize our notation, which follows that of Ref. [JX3 



and we introduce and describe key properties of several superoperators that are important 
for our analysis. 

The space of linear operators acting on a Hilbert space 7i is a Z) 2 -dimensional complex 
vector space. We introduce operator "kets" \A) = A and "bras" (A\ = A\ distinguished from 
vector kets and bras by the use of smooth brackets. The natural operator inner product can 
be written as (A\B) = ti(A^B). An orthonormal basis je^) induces an orthonormal operator 
basis 

\ej)(e k \ = r jk = r a , (Al) 

where the Greek index is an abbreviation for two Roman indices. Not all orthonormal oper- 
ator bases are of this outer-product form. In the following, r a can be a general orthonormal 
operator basis, or it can be specialized to an outer-product basis. 

The space of superoperators on TC, i.e., linear maps on operators, is a Z) 4 -dimensional 
complex vector space. A superoperator A is specified by its "matrix elements" 



A, 



lj,mk 



*4(|ei>(e fc |) e m ) , (A2) 
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for the superoperator can be written in terms of its matrix elements as 



A= E •Aij,mk\ei){e j \ \e k )(e m \ = E t q = A a/ 3\r a ) (r^ . (A3) 

lj,mk a,/3 a,/3 

The ordinary action of A on an operator A, used above to generate the matrix elements, is 
obtained by dropping an operator A into the center of the representation of A, in place of 
the sign, i.e., 



A{A) = J2 A a( s t*AtI . 

There is clearly another way that A can act on A, the left-right action, 

A\A) =J2A a f3\T a )(Tp\A) , 



in terms of which the matrix elements are 



A a p = (r Q | A\rp) = (|ej)(ej 



A 



ei 



A(\ej)(e k \ 



&m / A, 



lj,mk 



(A4) 



(A5) 



(A6) 



This expression provides the fundamental connection between the two actions of a superop- 
erator. 

With respect to the left-right action, a superoperator works just like an operator. Mul- 
tiplication of superoperators B and A is given by 



BA= E Ba-yA 7 p\T a )(Tp\ , 

and the "left-right" adjoint, defined by 

(A\A\B) = (B\A\A)* , 

is given by 

a ] = T l = E-^MM • 

a,f3 a,f3 



(A7) 



(A8) 



(A9) 



With respect to the ordinary action, superoperator multiplication, denoted as a composition 
B o A, is given by 



B o A = E B-isAxp T ~/ T a t\tI . 

The adjoint with respect to the ordinary action, denoted by A x , is defined by 

tr([A x (B)]^A) =tr(5U(A)) . 
In terms of a representation in an operator basis, this "cross" adjoint becomes 

A X =E^/3 r l© T /3- 

a,/3 



(A10) 



(All) 



(A12) 
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Notice that 



(B o A) 1 = B ] o A ] and (BA) x =B X A X . (A13) 

We can formalize the connection between the two kinds of action by defining an operation, 
called "sharp," which exchanges the two: 

A*\A)=A(A). (A14) 

Simple consequences of the definition are that 

(A*^ = {A x 'f , (A15) 
(Bo A)* = B*A* . (A16) 

The matrix elements of A* are given by 

*4*mfc = (|e«)(e j ||^l # ||e m )(e fc | 



tr 



(| ej )(e^(|e m )(e fc |)) 
-A(|e m )(e fe |) < 



ei 



which implies that 



Alm,jk , (A17) 



A* = Aj,mk\ei)(e m \ |e fc > <e J -j . (A18) 

lj,mk 



A superoperator is left-right Hermitian, i.e., A^ = A, if and only if it has an eigende- 
composition 

A = J2fi a \r a )(r a \ = ^2^ a T a QTl, (A19) 

a a 

where the real (left-right) eigenvalues and the operators orthonormal eigenop- 

erators. 

A superoperator is trace preserving if, under the ordinary action, it leaves the trace 
unchanged, i.e., if tr(A) = tr (A(A)) = tr([A x (I)y A) for all operators A. Thus A is trace 
preserving if and only ii A x (I) = I . 

A superoperator is said to be positive if it maps positive operators to positive operators 
under the ordinary action. A superoperator is completely positive if it and all its extensions 
X ® A to tensor-product spaces, where X is the unit superoperator on the appended space, 
are positive. It can be shown that A is completely positive if and only if it is positive 
relative to the left-right action, i.e., (A|^4|A) > for all operators A (for a proof in the 
present notation, see Ref. flI5| ). This is equivalent to saying that A is left-right Hermitian 
with nonnegative left-right eigenvalues. 

In this paper we make use of several special superoperators, whose properties we sum- 
marize here. The identity superoperator with respect to the ordinary action is 

X = /0/ = ^|e,)(e j |0|e fc )(e fc | . (A20) 
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This superoperator is Hermitian in both senses, i.e., X = X' = X x . It is the identity 
superoperator relative to the ordinary action because 1{A) = A for all operators A, but its 
left-right action gives T\A) = tr(A)I. 

The identity superoperator with respect to the left-right action is 

1 = S K)(r a \ = \ej)(e k \ © |e fc )<ej| . (A21) 

a j,k 

This superoperator is also Hermitian in both senses, i.e., I = 1^ = I x . It is the identity 
superoperator relative to the left-right action because 1\A) = A for all operators A, but its 
ordinary action gives 1(A) = tr(A)I. Since sharping exchanges the two kinds of action, it is 
clear that X* = I. 

To define the remaining superoperators, it is useful to introduce a set of D 2 — 1 tracefree, 
Hermitian operators |T7j], which are the generators of SU(D). We label these operators by 
a Greek index a, which runs from 1 to D 2 — 1. The operators are defined by 

a = 1,...,D- 1 : 

Aa = Ij= , = ( Er kk -(j-l)r n ) , 2<j<D, (A22) 
~ 1) \k=i ) 

a = D,...,(D + 2)(D-l)/2: 

K = = 4| (Tifc + ^i) , 1 < i < k < D , (A23) 
a = £)(£) + 1)/2,...,D 2 - 1 : 



A Q = ry = -j= ( Tjk - r kj ) , 1 < j < k < D . (A24) 

In Eq. ( |A22| ), a stands for a single Roman index j, whereas in Eqs. (|A23| ) and (|A24 ), it 
stands for the pair of Roman indices, jk. These operators are Hermitian generalizations of 
the two-dimensional Pauli operators: the operators ( |A22| ) are diagonal in the chosen basis, 
like a z ; for each pair of dimensions, the operators (|A23|) are like the Pauli operator a x ; and 
for each pair of dimensions, the operators (A24) are like a y . 

Like the Pauli operators, the operators A Q are orthonormal, i.e., 

(X a \Xf3) = tr(A a A^) = 5 a p . (A25) 

Thus they constitute an operator basis for the subspace of tracefree operators. Indeed, we 
can define a superoperator projector, 

^ee^|A q )(A q | = £A Q OA Q , (A26) 

a a 

which relative to the left-right action, projects onto the subspace of tracefree operators. 
Notice that T = P = J rx . 

If we add to the set of operators X a the normalized unit operator I/yD, we obtain an 
orthonormal operator basis. Thus the unit superoperator I can be written as 

I = + £ |A Q )(A Q | = 1/D + T . (A27) 



—i 
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Writing T = I — 1/ D, we find that 



I 

D 



D 



— 

£> 2 X> 



(A28) 



In the chosen basis, the operators ( |A22| ) and (A23) are real and symmetric. Together with 
I/y/D, they constitute a set of D(D + l)/2 orthonormal operators, which span the subspace 
of operators that are symmetric in the chosen basis. In contrast, the D(D — l)/2 operators 
in Eq. (|A24| ) are pure imaginary and antisymmetric and span the subspace of operators that 
are antisymmetric in the chosen basis. We can define superoperator projectors, 



'Ps^^P+ £ | Ac) (Ac 



D 



A a imaginary 



A a real 

|A a )(A c 



(A29) 
(A30) 



which relative to the left-right action, project onto the symmetric and antisymmetric oper- 
as — >'S — r S miu r A — i' A — r A • 



ator subspaces. Notice that Vs = V\ = Vg and Va = V\= Va ■ It is clear that 



l = V, 



s 



(A31) 



The last superoperator we need is the superoperator that transposes operators in the 
chosen basis. The ordinary action of the transposition superoperator is given by 



T{A) = Y J \e J ){e k \A\e J )(e k \ 



so the superoperator has the form 



T = E l e i)( e fcl © l e i)( e fcl • 



(A32) 



(A33) 



The transposition superoperator is Hermitian in both senses and is unchanged by sharp- 
ing, i.e., T = T' = T x = T*. In addition to satisfying T o T = X, the transposition 
superoperator has the property that 



I o T = I , (A34) 

which in view of Eq. ( |A16| ) , is equivalent to IT = I. 

It is easy to see that Vs — Va, acting to the right, transposes an operator, i.e., 

V S \A) - V A \A) = T{A) = T*\A) , (A35) 

which gives us, since T is invariant under sharping, 

T = T* = V S -V A ■ (A36) 

Combined with Eq. (|A31|) , this gives us 
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V S = -(I + T), 
V A = 1 -(I-T). 
Combining these forms with Eq. (|A34|) yields 



2V s oT = 1 + 1 = (D + l)g AV 
2V A o T = I - J = S D /v D . 



(A37) 
(A38) 



(A39) 
(A40) 



The form ( A40 ) has been given previously 



APPENDIX B: 



In this Appendix we show that a superoperator is Hermitian relative to the left-right 
action if and only if it maps all Hermitian operators to Hermitian operators. 

Let A be a superoperator, and let \ej) be an orthonormal basis, which induces an or- 
thonormal operator basis |ej)(efc|. Notice that 



ei 



AH\e,)(e k \) 



( |e m )(e fc | 



A 



ei 



A\e k )( ej \) 
iA\e k )(e 3 \)V 



ei)(ej\ 



(Bl) 



Here the first and third equalities follow from relating the ordinary action of a superoperator 
to its left-right action [Eq. (|A6|) 1, the second equality follows from the definition of the left- 
right adjoint of A [Eq. (|A8|)], and the fourth equality follows from the definition of the 
operator adjoint. Equation ( [Bl|) gives the relation between the operator adjoint and the 
left-right superoperator adjoint: 



^(| e .)( efc |) = [^(| efe )( e ,|)]t. 
Thus we have that A = A\ i.e., A is left-right Hermitian, if and only if 



^(|e J -><c*|) = [^(|e*>< 



It 



(B2) 



(B3) 



for all j and k. This result allows us to prove the desired theorem easily. 

Theorem. A superoperator A is left-right Hermitian if and only if it maps all Hermitian 
operators to Hermitian operators. 

Proof: First suppose A is left-right Hermitian, i.e., A = A^ . This implies that A 
has a complete, orthonormal set of eigenoperators r a , with real eigenvalues fi a . Using the 



eigendecomposition ( |A19| ), we have for any Hermitian operator H, 



A(H) = J2^r a Hri = A(Hy. 



(B4) 
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Now suppose A maps all Hermitian operators to Hermitian operators. Letting Tj 



|ej)(efc|, it follows that 



A{r jh 



«M ^( T jk + r kj ) + i-^(r jk - T kj 



A[ \i r jk + T kj )) - iA(-^(r jk - r kj ) 

it 



A(^-{T jk + T kj ) - i — (T jk ~ T kj ] 

[A(r kj ^ . 



(B5) 



Equation (|B3|) then implies that A — A'. ■ 
Since a superoperator is left-right Hermitian if and only if it has an eigendecomposition as 
in Eq. (|A19| ), we can conclude, by grouping together positive and negative eigenvalues, that 
being left-right Hermitian is equivalent to being the difference between two completely pos- 
itive superoperators. Using the theorem, we have that a superoperator takes all Hermitian 
operators to Hermitian operators if and only if it is the difference between two completely 
positive superoperators. This generalizes a result of Yu fI5|| , who showed that a positive 



superoperator is the difference between two completely positive superoperators. From our 
perspective, we can say that since a positive superoperator takes positive operators to posi- 
tive operators, it also takes Hermitian operators to Hermitian operators and thus is left-right 
Hermitian. A positive operator that is not completely positive has one or more negative 
left-right eigenvalues. 

We can get one further result relevant to the considerations in this paper: if A and B are 
left-right Hermitian superoperators for two separate quantum systems, then A ® B is also 
left-right Hermitian and thus maps all Hermitian operators of the joint system to Hermitian 
operators. 



Let 
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D 1 D 2 

PA = J2^j\ej)(ej\ and p B = Vk\fk)(fk\ 

j=l k=l 



(CI) 



be the eigendecompositions of pA and ps- In the joint basis \e 3 ) ® \ f k ), pab has the form 

PAB = Pjk,lm\ e j)( e l\®\fk)(fm\- (C2) 

j,k,l,m 



The diagonal forms of the marginal density operators show that 
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D 2 Di 

^ Pjk,lk = P>j$jl and X! Pjkjm = ^fc^fcm • (C3) 

fe=l i=l 

Thus the diagonal elements of Pjk,i m are a probability distribution pj k = Pjkjk, whose 
marginals are the eigenvalues of the marginal density operators: 

D 2 Di 

^Pjk = Pj and J2pjk = Vk- (C4) 

fc=l j=l 

We now can write 



1 +^(Pab) = 1 + 53 \Pjk,lm\ 
j,k,l,m 



2 



53 PjkPlm + J2P% 
j,k,l,m j,k 

53 + 53 PjfePim + 53 pjkPik - 53 pjkPik 

j,k,m jjtl,k,m j,k,l jj^^k 

Z)fe Pi*) + 53 (53 Pi*) + 53 Pi*Pftn 

i fc fc i j=£l,k=£m 



>£p? + 53 



ft 



i i 

= tr(p^)+tr(p|). (C5) 

The first inequality here is saturated if and only if p AB is diagonal in the basis \e 3 ) ®\fk)- 
The second inequality is saturated if and only if PjkPim = whenever j 7^ I and k 7^ m. 
This requirement is equivalent to saying that the nonzero entries in pjk are restricted to one 
row or to one column. In view of the first requirement, this means that overall equality is 
achieved in Eq. (|U5|) if and only if p AB = p A <g> p B is a product state, with p A or p B a pure 
state. 



APPENDIX D: 

In this Appendix we show that the vector space of operators acting on a D-dimensional 
Hilbert space has only two proper operator subspaces that are invariant under all unitary 
transformations. These two subspaces are the one- dimensional subspace spanned by the 
unit operator / and the subspace consisting of all tracefree operators. 

It is obvious that the subspace consisting of multiples of / and the subspace of trace-free 
operators are unitarily invariant. To show that these are the only unitarily invariant proper 
subspaces, we consider a unitarily invariant subspace that is not the subspace spanned by 
J, and we show that this subspace is either the subspace of tracefree operators or the entire 
operator space. Let A be a nonzero operator in the unitarily invariant subspace, which is 
not a multiple of /. There exists an orthonormal basis \ej) such that An 7^ ^22- Adopt this 
basis, in which A has the representation 
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A=J2 A jk \e s )(e k \. (Dl) 
j,k=l 

Consider the unitary operator U that changes the sign of |e x ) , i.e., U\e-\) = — |ei) and 
U\ej) = \ej) for j = 2, . . . , D. Also in the unitarily invariant subspace is the operator 

1 D 
B = -(A+UArf)=A 11 \e 1 )(e 1 \ + £ A jfc | ej -)(e fc | . (D2) 

^ j,fc=2 

Do the same thing to the second basis vector; i.e., use the unitary operator V defined by 
V\ e 2) = — |e 2 ) , and V\ej) = \ej) for j = 1 and j = 3, . . . ,D. Also in the subspace is the 
operator 

1 D 

C = -(B + VBV^)=A 11 \e 1 )(e 1 \+A 22 \e 2 )(e 2 \ + £ A,- fc | ej )<e fc | . (D3) 

Z j,k=3 

Now consider the unitary operator W that swaps \ei) and |e2), i.e., VF|ei) = \e 2 ), W\e 2 ) = 
|ei), and W|ej) = \ej) for j = ?>,..., D. Also in the subspace is the (nonzero) tracefree 
operator 

D = C — WCW^ = (A n - A 22 )(| ei )( ei | - |e 2 )(e 2 |) . (D4) 

We conclude that the subspace contains the tracefree operator |e 1 )(e 1 | — |e 2 )(e 2 |, which 
is a Pauli a z operator for the first two dimensions. From this operator, we can generate by 
unitary transformations that interchange basis vectors a <j z -like operator for every pair of 
dimensions, and from these a z operators, we can generate by unitary transformations a a x 
and a a y operator for every pair of dimensions. Since these Pauli-like operators span the 
space of tracefree operators, we conclude that any unitarily invariant operator subspace that 
is not the space spanned by / contains all tracefree operators. 

The unitarily invariant subspace could be the subspace of tracefree operators. Suppose 
that it is not and thus contains an operator E that is not tracefree. Defining a tracefree 
operator F = E — tr(E)I / D, we see that I can be written as linear combination of F and 
E and thus is in the subspace. Since the tracefree operators together with / span the entire 
space of operators, we conclude that in this case the unitarily invariant subspace is the entire 
operator space. This establishes our result. 
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